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By considering correlations between classical orbits we derive semiclassical expressions for the 
decay of the quantum fidelity amplitude for classically chaotic quantum systems, as well as for its 
squared modulus, the fidelity or Loschmidt echo. Our semiclassical results for the fidelity amplitude 
agree with random matrix theory (RMT) and supersymmetry predictions in the universal Fermi 
golden rule regime. The calculated quantum corrections can be viewed as arising from a static ran- 
dom perturbation acting on nearly self-retracing interfering paths, and hence will be suppressed for 
time-varying perturbations. Moreover, using trajectory-based methods we show a relation, recently 
obtained in RMT, between the fidelity amplitude and the cross-form factor for parametric level cor- 
relations. Beyond RMT, we compute Ehrenfest-time effects on the fidelity amplitude. Furthermore 
our semiclassical approach allows for a unified treatment of the fidelity, both in the Fermi golden 
rule and Lyapunov regimes, demonstrating that quantum corrections are suppressed in the latter. 

PACS numbers: 03.65.Sq, 05.45.Mt 



I. INTRODUCTION 

Understanding the sensitivity of quantum time evolu- 
tion in complex systems with respect to perturbations 
has evolved to a common subject of interest in the com- 
plementary fields of quantum information and quantum 
chaos. The concept of fidelity has become a central mea- 
sure to quantify the stability of quantum dynamics upon 
perturbations. The fidelity M{t) was introduced by Peres 
[l| as the squared modulus of the fidelity amplitude, 
m(t), the overlap integral of an initial state, e.g. a wave 
packet, with the state obtained upon forward and back- 
ward propagation governed by two Hamiltonians differ- 
ing slightly by a perturbation. In the context of quantum 
chaos the fidelity is often referred to as Loschmidt echo 
2] , picking up and generalizing concepts from spin echo 
physics Q. 

By definition the fidelity M(t) equals unity at t = 
and usually decays further in time. Broad interest in 
fidelity decay had been initiated by the pionieering semi- 
classical work of Jalabert and Pastawsky 0] who found a 
perturbation-dependent decay for weaker perturbations 
and discovered the intriguing Lyapnov regime for clas- 
sically chaotic quantum systems under stronger pertur- 
bations. Nowadays, three main decay regimes are dis- 
tinguished, depending on the strength of the perturba- 
tion S: In the limit of a weak perturbation, i.e. if £ 
is smaller than the quantum mean level spacing A, the 
fidelity decay is Gaussian characterizing the perturba- 
tive regime [y, @, 0|- For perturbations of the order of 
or larger than A, the decay is predominantly exponen- 
tial, M(t) ~ e~ 2rt , with a decay constant T obtained 
by Fermi's golden rule. The corresponding perturba- 
tion range is hence called the Fermi golden rule (FGR) 
regime 0, H, @. For strong perturbations the decay is 
still exponential in time, but with a ^-independent decay 



rate given by the Lyapunov exponent A of the classically 
chaotic counterpart of the quantum system characteriz- 
ing the Lyapunov regime 0, Q ■ For a review of the exten- 
sive literature on fidelity decay see Refs. [UEJ], including 
a discussion of further decay mechanisms. 

Fidelity decay has been calculated within two main 
theoretical frameworks, namely random matrix theory 
(RMT) and supersymmetrical approaches on one side, 
and semiclassical theory on the other. Within RMT, the 
decay of the averaged fidelity amplitude was computed in 
the perturbative and FGR regime within linear response 
ill] for weak perturbations. Later, using supersymmet- 
ric techniques and going beyond linear response, these 
results were extended to strong er p erturbations for the 
different universality classes [12|, LL3| • For the GUE case, 
the result for the fidelity amplitude in the universal FGR 
regime takes the particularly compact form 
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with 7 = Tt, t = t/Tn and T H = 2nhA is the Heisen- 
berg time. In particular, these supersymmetrical cal- 
culations revealed a fidelity recovery at the Heisenberg 
time. Recently, this has been interpreted by establish- 
ing relations between the fidelity amplitude decay and 
parametric level correlations [l4| . allowing fidelity decay 
phenomena to be considered from a completely different 
perspective. There it was found that in the FGR regime 
the Fourier transform of the parametric level correlator, 
the cross-form factor K (r;7), and the RMT fidelity am- 
plitude to(t; 7) are closely related to each other through 
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for the different GOE. GUE and GSE RMT ensembles 
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labeled by the index (3—1,2 and 4. 

Despite this recent progress, the averaged fidelity M(t) 
itself, the Loschmidt echo, has not been obtained using 
supersymmetric techniques. Within an RMT framework, 
M(t) was calculated within linear response which - upon 
heuristically exponentiating this result - gives a fair ap- 
proximation for the transition region between the per- 
turbative and FGR regime [HI- Both, RMT and su- 
persymmetry approaches are limited in principle to the 
perturbative and universal (FGR) regime, governed by 
a single, system-independent perturbation parameter F, 
and cannot reveal information about the individual sys- 
tem, such as the Lyapunov decay for strong perturba- 
tions. Hence, the Lyapunov regime is not amenable to 
RMT approaches. 

The complementary situation appears from the view- 
point of semiclassical theory 0, H, [3 EBl : While the ex- 
isting semiclassical tools have lead to the discovery of the 
Lyapunov regime, M(t) ~ e _At [H, uncovering an appeal- 
ing connection between classical and quantum chaotic 
dynamics, semiclassics has so far only been able to pre- 
dict the leading exponential decay in the FGR regime [4| , 
but could not account for additional quantum corrections 
arising from an expansion, in orders of t/Tn, of the super- 
symmetry results such as Eq. (TTJ). The absence of such 
quantum interference contributions in the semiclassical 
approach can be traced back to the so-called diagonal 
approximation used in the treatment of the fidelity am- 
plitude m(t), that is a pairing of the same trajectories in 
a semiclassical path integral representation of the prop- 
agators for the perturbed and unperturbed system. 

In this paper we will go beyond this approximation by 
evaluating contributions from correlated trajectory pairs 
built from different orbits. Classical correlations between 
(periodic) orbits were shown to be the key to understand- 
ing RMT predictions for spectral statistics [l7| and to 
deduce universal spectral properties within a semiclas- 
sical theory. This technique has been considerably fur- 
ther develo ped and applied to calculate various spectral 
0, IS |2(1 21 22 , 23 1 ,scattering and transport proper- 
ties [2J, 12a. 2a , [271 |28[ . [29j of classically chaotic quantum 
systems [30 ■ Recently, this method has been extended 
to semiclassically compute the quantum survival proba- 
bility and photo-fragmentation processes of open chaotic 
systems [111, HH, as well as to establish a semiclassical 
version of the continuity equation [33j . Classical correla- 
tions encoded in the semiclassical diagrams considered in 
the works [3l|, HU, [H| are of special relevance for ensuring 
unitarity in problems involving semiclassical propagation 
along open trajectories inside a system and will prove 
particularly important for the fidelity decay. 

We will show in SeclUfor the fidelity amplitude in the 
FGR regime for times below Th that such subtle classi- 
cal correlations also provide the semiclassical key to the 
aforementioned quantum corrections to the exponential 
decay of m(t), arising from an expansion of Eq. (jTJ) in 
powers of t/T-R. It is moreover of interest to explore the 
implications of the Ward identities leading to Eq. @ on 



the level of the semiclassical theory. In Sec. IIIII we de- 
rive this relation by invoking recursion relations for terms 
containing the correlations between classical orbits. 

Beside the formal derivation of Eq. ([TJ) , semiclassics 
also sheds light on the underlying interference mechanism 
leading, e.g., to weak- localization- type corrections to the 
exponential fidelity amplitude decay, encoded in Eq. (jTJ) . 
This arises from the fact that the relevant dynamics is or- 
ganized along orbits with so-called encounter stretches, 
orbit segments which are traversed (time-shifted) twice in 
nearly opposite directions. On its way back and forth the 
particle experiences practically the same disorder pertur- 
bation potential. This gives rise to an enhanced single- 
particle damping compared to dynamics along a generic 
orbit exploring independent disorder regions, and hence 
finally leads to an enhanced fidelity decay. This has inter- 
esting physical implications which we will discuss below: 
(i) such interference effects will depend on the Ehrcnfcst 
time where the correlation length of the disorder poten- 
tial appears as a new relevant length scale as will be 
discussed in Sec. IIVI (ii) The above outlined investiga- 
tion refers to a static (random) perturbation as most of 
the works on fidelity decay. However, deviations are ex- 
pected for time-varying perturbations, as e.g. caused by 
phonons. Our approach allows for incorporating such 
effects. We will show in Sec. [V] how a decrease in the 
correlation time, i.e. faster time variations, gives rise to 
a stronger suppression of the quantum interference terms 
and thereby in turn to an increase of the overall fidelity 
amplitude. 

In Sec. |Vl]we address the fidelity M(t), which is usu- 
ally easier to measure. We compute the corresponding 
quantum corrections for M(t) and show that they are 
suppressed in the Lyapunov regime. Our semiclassical 
approach hence adequately describes the FGR and Lya- 
punov fidelity regimes in a unified way. 



II. FIDELITY AMPLITUDE 
A. Semiclassical approach 

The fidelity amplitude m(t) is defined as the overlap 



m{t) 



dvS>* (r, i)$(r,i), 



(3) 



between states "J/ (r, t) and <!> (r, t) obtained by propagat- 
ing a common initial state \to with two slightly differ- 
ent Hamiltonians: H and H' = H + Yj. The integration 
in Eq. ([3]) is performed over the available configuration 
space, which in the case of billiards is the area enclosed 
by the boundaries. For further semiclassical analysis of 
m(t) it is convenient to represent ^ (r, t) with the help of 
the quantum propagator Kjj (r, r', t) for the Hamiltonian 
H: 



tf (r,t) = J dr'K H (r,v',t) * (r') 



(4) 
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The same representation holds for <f>(r, i), with Kjj re- 
placed by the propagator K^i for the perturbed Hamil- 
tonian H' . The next step is to use the semiclassical Van 
Vleck formula [35[ for the propagators K h and Kjj 1 ■ For 
two-dimensional systems one has 



^sc(r,r',t) 



2-Kih 
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7(r'— *r,t) 



C-' exp 



-^(r,r',t) 



(5) 

where S^(r,r',t) is the classical action of the trajec- 
tory 7 running from r' to r in time t, and C 7 = 
|det (d 2 S^/dr'dr) | exp (—^^7) is the Van Vleck deter- 
minant including the Maslov index /i 7 . 

We will assume that the perturbation £ is classi- 
cally small such that only the actions, i.e. the phases, 
are affected while the classical trajectories 7 remain un- 
changed. Under this assumption, after inserting Eqs. 
and (|4]) into Eq. ([3]) we obtain the following semiclassical 
approximation [4j for m(t): 
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exp 
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where AS^ stands for the change in the action along 
the trajectory 7 due to the perturbation S. Note that 
Eq. ([6]) involves a double sum over trajectories 7 and 
7' of the unperturbed Hamiltonian H only. Due to the 
rapidly oscillating phase factor containing the action dif- 
ferences Sj — Sy' most of the contributions will cancel out 
except for semiclassically small action differences origi- 
nating from pairs of trajectories which are close to each 
other in the configuration space. We thus can use a lin- 
ear approximation in order to relate the actions 
along the trajectories 7, 7' to the actions S* 7 , S~ ( > along 
the nearby trajectories 7, 7' connecting the midpoint, 
r = (r' + r") /2 with r. This yields 



m sc (t) 



JL) /drdro«(r + |)*o(ro-f 



x exp 



AS- 



(7) 



where q = (r' — r"), and Pq and Pq are the initial mo- 
menta of the trajectories 7 and 7', originating from the 
expansion of the actions around ro . 



B. Treatment of the perturbation 

Further evaluation of m sc (t) requires us to consider the 
effect of the perturbation on the phase difference AS 7 in 



Eq. ([7]). It can be expressed as 



(8) 



where (if) is the difference between the kinetic and 
potential energy of the perturbation E. In the case of a 
perturbation potential S, Eq. J5J) simplifies to 



AS 7 (<) = - / di'£[q 7 (*'),*']■ 



(9) 



For fully chaotic systems the AS 1 (t)/t are distributed 
as Gaussian random variables when t is sufficiently large 
(compared to all classical time scales), and the variance 
is given by t d T C(r), where C(r) = (L= (r) (0)) 
is time correlation of the perturbation. Assuming that 
the mean value of L 7 is zero this implies 



exp ( ^AS" 7 



exp 



exp 



2h 2 



where 



(10) 



(11) 



Note that the average in Eq. (fTU|) should be understood 
as the average in the phase space over trajectories 7 with 
different initial conditions. Alternatively, one can fix 7 
and consider an average over an ensemble of different 
perturbations of L^. One possible realization of such 
an ensemble is provided by a quenched static disorder 
potential 



s w-E^ ex p 



(12) 



It is given by Ni random impurities in a cavity of an area 
A with the Gaussian profile characterized by a finite cor- 
relation length £, see [Ullij]. The independent impurities 
can be assumed to be uniformly distributed at positions 
R a with the densities n,; — Ni/ A and strengths u a obey- 
ing (UaUff) = u 2 S al3 . 

The action differences AS 7 accumulated by segments 
of 7 separated by distances larger than £ can be regarded 
as uncorrelated. Consequently the stochastic accumula- 
tion of AS-y along 7 can be described by a random pro- 
cess, resulting in a Gaussian distribution of the action 
difference AS" 7 . This yields the previous result (TIT))) with 
the noticeable difference that the avera ge h ere is over the 
ensemble of the disorder potentials [4|, [3J| . For the po- 
tential in Eq. (Ti2"]) the decay time t in Eq. (fTT)|) can be 
calculated explicitly. In this case the correlation function 
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is given by 



Csflr-r']) = (S(r(t))E(r(f'))> 
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(13) 



and it depends only on the difference between r = r(t) 
and r' = r(t'). Using then ergodicity of the classical flow 
and substituting time averages in Eq. (ITTj) by the inte- 
gration over the cavity domain we obtain for the decay 
time the relation I34| 



f£ = V2tt k£ t 5 



(14) 



Here, k = mvo/h = 2tt/Xb, with vq the particle velocity 
and Xb the de Broglie wave length, and is = h 3 /(niU 2 m) 
the quantum elastic scattering time for the white-noise 
case of (5-scatterers. Note that for k^ ^> 1, coincides 
with the clastic scattering time obtained quantum me- 
chanically within the first Born approximation for the 
disorder potential (jT2J) . In the limit fc£ < 1, where the 
semiclassical treatment of disorder effects in terms of un- 
perturbed trajectories is no longer valid, Eq. (JTUJ) can still 
be used but with replaced by tg. 



C. Evaluation within diagonal approximation 

In order to evaluate the double sum over paths in Eq. 
(|7|) let us first consider only the diagonal contribution by 
pairing each path with itself, i.e. arising from the terms 
7'. After inserting Eq. (|10p into Eq. ([7]) we obtain 



m d {t) 



= f JL 



\2irh 

\ C t\ e ^P 



r/r-Zi-u-Zq*, 1 , (r + |) #0 (r - | 



7(ro— >r,t) 



:QP6 exp 
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The prefactor |C 7 | can be regarded as a Jacobian when 
transforming the integral over the final position r into 
an integral over the initial momentum po = Pq of each 
trajectory. We then obtain 



m d (t) 



(16) 



Here 



ro po indicates the phase space average |32| . 
( f )r , P0 = ( 27 r/t)2 / dr odPoF(r ,p )pw (r ,Po), 



(17) 



where 



p w (r , po) = J dq% (r + |) *o (ro - | 



x exp [ --^qPo 
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FIG. 1: Sketch of pairs of correlated classical trajectories 7 
and 7'. In (a) the encounter region (rectangular box) connects 
a loop with two long legs beginning at ro and ending at r. In 
contrast the paths begin or end inside the encounter region 
('one-leg- loops') in (b). The zoom into the encounter region 
in (b) illustrates the position of the Poincare surface of section 
(PSS) used and the definition of the encounter time. 



denotes the Wigner function of the initial wave packet 
^o- Assuming that V^o has a small energy dispersion 
around a mean energy Eq, Eq. (|16p can simply be re- 
placed by 



m d (t) = exp (-It) 



(19) 



with the decay rate T = l/(2i (Eq)). This represents the 
exponential FGR decay in the universal regime 0, S @] • 
As our expressions for the fidelity amplitude depend on 
r, we will also use in the following the notation m (t; T) 
instead of m(t). 



D. Loop corrections 

We now consider contributions to Eq. |7]) from pairs of 
different trajectories following each other closely in con- 
figuration space, shown as full and dashed line in Fig. 
[T] The structure of such pairs can be characterized by 
long links, where two trajectories almost coincide, and 
by encounter regions (the boxes in Fig. [T]) where the seg- 
ments of two trajectories are connected in a different way. 
Correlated trajectory pairs of this type were first intro- 
duced in the context of spectral statistics involving peri- 
odic orbits by Sieber and Richter [17] ] showing how uni- 
versal RMT predictions can be obtained based only on 
semiclassical methods. We already mentioned important 
extensions of this approach in the introduction. Within 
the phase space framework developed in Refs. [l8| we will 
closely follow the lines of Refs. [31, 32], where this formal- 
ism was employed in a semiclassical approach to decay 
and fragmentation processes. The main difference to the 
present case consists in the fact that the quantum cor- 
rection to the classical exponential survival probability in 
Refs. [33, HH (due to trajectory pairs containing a com- 
mon encounter region) has to be replaced by a different 
correction to Eq. (fTB|) which has to be computed. One has 
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to account for the fact that the disorder induced phase 
difference is modified for correlated trajectory pairs, be- 
cause the same perturbation acts twice within the en- 
counter region. Although an expression for such a phase 
difference is known in a general context [2f|, we first 
illustrate how it can be calculated in the case of an orbit 
with the simplest possible encounter shown in Fig. [T^. 
To this end we split the orbit into parts consisting of two 
encounter-stretches and three links, 



exp 




(20) 



In the last line we assume that, in the semiclassical limit 
and for a fixed disorder correlation length £, the nearly 
parallel encounter stretches in the box are so close to each 
other (at a distance smaller than £) that they experience 
the same perturbation. Considering a Gaussian phase 
distribution for each orbit segment (encounter stretch 
and three links) and neglecting further correlations be- 
tween the different segments we finally apply Eq. (|f Op to 
every part of the orbit individually and obtain 



exp ( — A/S-y 



exp 



2h 2 

exp [-r (t + 2i cnc )] 



A<? 2 
A ^(i) 

2ft 2 



(21) 



Here t enc denotes the length of the encounter region. 

Apart from the basic diagrams in Fig. QJi we will also 
consider trajectory pairs that differ in an arbitrary num- 
ber of encounters with an arbitrary number of stretches 
involved. For quantifying the encounter structure of an 
orbit, we introduce notation as in Refs. [H, |26j. Wc 
define a vector v = (v2,V3, ...), where the component 
vi gives the number of Z-encounters, i.e. the number of 
times where I stretches of an orbit come close to each 
other. We denote by V = v i tne overall number of 

encounters and by L = ^4=2 the overall number of 
encounter stretches of an orbit. Repeating now the steps 
in Eqs. (|2"01|2"T1) for the general case of multiple encounters 
one finds 21. 12211 
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exp 
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a = l 



(22) 



where l a denotes the number of stretches of the encounter 
a. We note that the structure of the correction to the 



survival probability in Refs. [311 132| is slightly different, 



exp 



-r d (t-£(i a -i)c 



(23) 



where r<j, the inverse dwell time 1/r^, takes the role 
of r. Using Eq. (|2"Tj) . we can now evaluate in Eq. ([7]) 
the off-diagonal contributions to the sums over trajec- 
tory pairs differing at encounters in the same way as in 
Refs. [3l|, HH, HII- Here we sketch some details of these 
calculations. In the region around each Z-encounter we 
consider a Poincare surface of section and measure the 
differences Sj and Uj between the piercing points of one 
of the trajectories along the stable and unstable mani- 
folds, respectively. In these coordinates the duration of 
the encounter a can be expressed as |18| 



t" 



•In 



max, «, 



max,- 



(24) 



and the action difference between two trajectories is given 
by 



v l a -l 

0=1 j=i 



(25) 



Similarly to Refs. [32|, [33| we have to distinguish three 
different cases when the density of encounters is consid- 
ered: (A) encounters inside the trajectory (Fig. [Th.), (B) 
encounters at the beginning or the end of the trajectory 
(Fig. [TJa) and (C) different encounters at the beginning 
and at the end of the trajectory. 

For trajectories of type (A) and duration t whose en- 
counters are characterized by v, the corresponding weight 
is given by 



w t .A (s,v) = N(v) 



^ L - y n;=i*e 



(26) 



Here, fi is the volume of the available phase space, and 
N (v) is the number of trajectory structures, i.e. the num- 
ber of topological realizations of orbits with the vector v, 
see Ref. |f8j . (Note that, in contrast to Refs. [12, [HI, we 
include N (v) into the weight function.) For a trajectory 
with the beginning or end point inside of an encounter, 
the weight function can be conveniently expressed as 
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Finally, for trajectories with both beginning and end 
point inside of two different encounters we have 



w t .c (s, u) 



t - v / t a 



L-2 



1,3=1 
0^/3 



(L-2)!^n:=i^nc 



(28) 



Here M a (v) is defined as the number of possible ways 
to cut links between the encounters a and f3 divided by 
L in all periodic orbit structures described by the vector 
v, see Ref. [H, HI]. 

We can now present general expressions for the loop 
contributions to the fidelity amplitude resulting from 
the three weight functions in Eqs. (|26M28p . In view of 
Eqs. (HU), (US]) and Eq. f2"6). we obtain in case (A) the 
following contribution to the fidelity amplitude from tra- 
jectories of time t with encounters structure characterized 
by v. 




m A (t;T,v) = N(v) 



dsdu 



x exp 



-r t + 



E (£ - la) cue 1 
a=l I . 



-e«° 



,(29) 



Applying then the rule [181 ] that, after expansion of 
Eq. (|29p in t" nc , the only non- vanishing contributions 
come from the terms independent of t" nc , we get 



m A (t;T,v) = 
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T H r 



L-V 
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I d 



wn - 5 
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(rt) L cxp (-rt) 



(30) 



where we used Th = 2nhA and the definition of compo- 
nents vi. In a similar way we can derive from Eq. (I27p 
the contribution for the case (B): 



m B (i;I» = 2 
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N(v)J2 
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(31) 
exp (-rt) 



(L-l)! 



with an Zi-encounter at the beginning or end of the tra- 
jectory. Finally we obtain from Eq. (f2"5| for the case (C) 



m c (t;T,v) 
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r H r 



L-V 
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2 y (32) 

" ; t L - 2 cxp (-rt) 



(L-2) 



with an ^-encounter and an Z 2 -encounter at the begin- 
ning and at the end of the trajectory, respectively. 

The entire contribution from orbit pairs to m (t; T) is 
obtained from Eqs. (|30j|32|) by summing over all possi- 
ble vectors v and finally adding up three contributions 



FIG. 2: Diagrams with two 2-encounters. In the unitary case 
only the diagram (c) contributes to the fidelity amplitude. 



together with the diagonal term (HHJ): 

m sm (t;T) = m d (t; T) + \m A (t; T, v) 

+m B (t;T,v)+m c (t;T,v) . (33) 

We now compare the contributions to the fidelity am- 
plitude resulting from our semiclassical expressions (|30I 
122) with the RMT results pj,[13- 

For systems with time-reversal symmetry the leading 
order correction is due to the 2-encounter diagrams in 
Fig. [TJ From these diagrams we obtain as one main result 



i qi 

GOE 



(*;r) 



-^-)«-Ni»i-r/) . cwi 



This leading quantum interference correction leads to 
a reduction of the fidelity amplitude. Equation (|34|) 
indicates that quantum deviations of the fidelity from 
pure exponential decay become relevant at time scales 
t ~ (iTn) 1 / 2 , with decay time t — \/2T, i.e. at times 
which depend on the perturbation strength and can be 
much shorter than the Heisenberg time Th- Thus the 
quantum corrections can occur at timescales before sat- 
uration sets in [J H, [3(| . 

The next-to-leading order contributions come from the 
diagrams of higher order with either two 2-encounters (as 
in Fig. [5]) or with one 3-encounter (as in Fig. [3]) [Hj]. All 
together they give: 



l GOE' 



(t; T) = ( -2—^ + ) exp (-Ft) . (35) 

H H 



In the unitary case the leading order correction origi- 
nates only from the diagrams shown in Fig. [2}; and Fig. 
|3K, yielding 



m q GUE (t; T) = ( -^y + ^ ) exp (-Ft) . (36) 



The results ()34ti36p are consistent with the ones obtained 
from RMT-calculations, as can be seen by expanding 



7 



<: ... 




FIG. 3: Diagrams with one 3-encounter. In the unitary case 
only the diagram (a) contributes to the fidelity amplitude. 



in t/Tptbe RMT expressions for m(t;T), t/T K < 1 in 
Refs. [H, [H and comparing the leading-order terms. 

It is instructive to analyze these results and the under- 
lying semiclassical assumptions for the case of the disor- 
der potential (fT2")) with correlation length £. To generate 
random (Gaussian) fluctuations, £ should be smaller than 
the system size C. Furthermore, in the semiclassical limit 
£ ^> Xb- During the derivation of the above equations we 
implicitly assumed that the distance between encounter 
stretches, of the order of ^/CXb, is smaller than £ in or- 
der to have correlated disorder along the stretches and 
thereby loop interference corrections. This implies that 
our approach is valid for parameters obeying 



V£A^<£ <£. 



(37) 



For white noise disorder, £ — > 0, the semiclassical ap- 
proach would predict vanishing interference corrections, 
since then the disorder potential along the two encounter 
stretches is no longer correlated. This would imply clear 
deviations from the universal RMT result for the white 
noise case. In that case the present semiclassical pic- 
ture of orbits which remain unaffected by the disorder 
(up to the phases) has to be replaced by a semiclassical 
approach based on paths scattered at impurities. We be- 
lieve that presumably Hikami boxes take the role of the 
encounter regions finally establishing universality for the 
white noise case. 



III. IDENTITIES FOR THE FIDELITY 
AMPLITUDE 



In this section we address two important properties 
of the semiclassical quantum fidelity amplitude m sc and 
prove them within the trajectory-based formalism for 
times t < Th- The first one: m sc (t;T — 0) = 1, is the 
direct consequence of the unitarity of the semiclassical 
evolution. The second one is the connection © between 
the fidelity amplitude and the parametric spectral form 
factor. As we will show below, from the semiclassical 
point of view both properties can be attributed to the 
existence of certain recursion relations satisfied by N (v) 
and Afki (v). 



The fact that we can prove to all orders both the uni- 
tarity of semiclassical evolution and Eq. ^ shows the 
consistency of our semiclassical approach for times below 



A. Unitarity 

In the case of vanishing perturbation (r = 0) the fi- 
delity amplitude m(t) should be equal to one by the uni- 
tarity of quantum evolution. In the following we will 
show that this property holds for the semiclassical form 
of m(t) obtained in the previous section when t < Th- 

As one can immediately see from Eq. (|T9")) . m sc (t; T — 
0) = 1 within the diagonal approximation. Hence one 
has to show that all further semiclassical loop contri- 
butions vanish for T = 0. To this end we will demon- 
strate that the off-diagonal terms niA(t; 0, v), msit; 0, v) 
and mc(t;0,v) cancel each other. According to 
Eqs. (|3QI31l32p they read: 



m A (t;0,v) = ^(-)(-l) y (4) L "-^^ ! ' 



m B (t;0,v) = 2N{v){-l) y - 1 ^-{l- 

x (L-vy. ' 

Y> (v) f J_ 

^ Ik 1 ' It h 

k,l=2 v 



(L-V)\ 



L-V 



(38) 



(39) 



m c (t;0,v) 



L-V 



(40) 



Adding these contributions, one can see that the off- 
diagonal corrections disappear if for each v the following 
condition is satisfied: 



N(v) - 2—N(v) + 
L 



y> A4i (v) 
^ kl 

k,l=2 



= 0. 



(41) 



Precisely the same unitarity condition was obtained in 
Ref. [1^ for the continuity equation. As it was shown 
there, Eq. (|4"Tj) can be proven by relating Mm (v) to N (v) 
via the relationship 

Nu[v) = (fc + ^Y L K _X 1 + 1) ^ W '' 11 )- ( 42 ) 

where Vk+i-i is the (fc + I — l)-th component of v, and 
v [k,i^k+i-i] j g .j-kg vec t or obtained by decreasing the com- 
ponents Vk and vi by one and increasing the component 
v k+i-i by one. To obtain Eq. (|42p one considers link con- 
traction such that a /c-encounter and a Z-encounter merge 
into a (k + l — l)-encounter. By looking at the number of 
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possible ways to contract the link and to form a smaller 
periodic orbit structure, Eg. can be deduced from 
the results in Refs. [11,01. 

In view of Eq. (|4"2)l . Eq. (|41[) can be transformed into 



L-2V-J2v k +i-i )N{v) =0 



(43) 



which can be proven after performing the double sums 
over k and I 13311. 



B. Relation between the fidelity amplitude and the 
parametric spectral form factor 

The surprising connection between the quantum 
fidelity amplitude and parametric spectral correlations 
was derived in Ref. [3| within an RMT approach. The 
idea was to use a certain invariance property of the inte- 
gration measure for the ensemble of random matrices. A 
corresponding Ward identity then led to Eq. (J2]), which 
can equivalently be put into the form 

T H /3 dK (t;T) 



= m (t; T) 



(44) 



2t 2 dT 

Here we show how this relationship can be obtained in the 
framework of our semiclassical approach for systems with 
and without time reversal symmetry ((3 = 1,2). To this 
end it is convenient to work with the Laplace transforms 
of the semiclassical expressions for the form factor and 
fidelity. 

In order to reveal a systematic structure for the con- 
tribution of each encounter stretch and each link, it is in- 
structive to take the Laplace transform of Eq. (f33|) with 
respect to 7 = Tt while keeping i] = tj (7h7) = (7H~r) -1 
fixed: 

F (<?: V) = J d>ym snl (jjThj; (^TnT 1 ) exp (-q-f) . 

Inserting the expressions for m sm (i; T) and performing 
the Laplace transformation by partial integration gives 
rise to the perturbative expression (in powers of nf) 



F(q,r ] ) = J2F n ( q )r 1 n , 



(46) 



where the n-th term in this expansion originates from 
trajectory pairs with L — V = n. Explicitly, the terms 
F n (q) take the form 

n£,H?-iy ~ 



F n (q) 



E 



L+l 



+2^2 N{v 



(9+1) 
l\Vh 



h=2 



N(v) 

g + 1 
L \-hq-ll 



E ,_/ , 



M lh («) (q+l) 



11) (■ 



-hq-ll) 



(47) 



where the sum is over all diagrams with fixed \ v\ = L — V. 
In Eq. (|47[) we recognize the following diagrammatic rule: 
every link contributes to F n (if) a factor (q+l) 1 and 
every Z-encounter a factor —l(q + I). 

Now we turn to the parametric spectral form factor 
for which we use the following semiclassical expression 
MM: 



K(t-Y) 



2Tt 2 ^f 1 \ L ~ V ( I d 



n 



/3T H ^V ri W 11 V tar 

v 1=2 



V 



xN(v) 



L — l -Tt 



IT 



(48) 



After taking the Laplace transform (|4"5|) of the left hand 
side of Eq. ([44]) we obtain for the n-th term of the ex- 
pansion (in powers of rj) 



d 



\v\-- 



N(v,q) 
L ' 



where N(v,q) is defined as 

N(v)(-l) v 



N(v,q) 



(49) 



(50) 



Upon performing the derivatives in Eq. (|49[) we obtain 



K= E 



N(v,q) 



\v\—n 



L 



qvi 

(q + l) 



(51) 



This must be compared with the F n term, Eq. (|47| . for 
the Laplace transform of the fidelity amplitude. Using 
Eq. ([50]) we can rewrite it as 




Furthermore, we can simplify this by expressing the ma- 
trix elements Nki (v) in terms of N(v, q). Using Eq. (j4"2| 
and taking into account the additional q-dependent fac- 
tors we get 



Fn = (-1) 



V-2 



Mu(v) Ui(i(q + i)Y 



k(q + k)l(q + l) (q + 1)^ 1 

N(v^«+^,q) vtir-1 +l ~ 1] 



(L-l) 



(q + k + l-iy 



(53) 



We can then rewrite the sum over the dummy vectors, 
v 1 = v l k , l ^ k + l ~ 1 ] ) ag a sum over v which gives for the 
second line in Eq. ([52]) : 



E 

\v\=r 



N(v,q) 
L 



E 



Vk+l-l 



(q + (k + l-l)) 



(54) 
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Using then the results of Ref. [33| and performing the 
sum over k it can be further simplified to 



V N(v,q) ( v (I - 2) Vl 
leading to the following expression for F n : 



(55) 



F„ 



E 

|v|=r 



N(v,q) 



(9 + 1) 



E 



Ivi 



(q + l) 



(56) 



The final step is to show that Eq. ([56]) coincides with 
Eq. ([51]) . This is indeed so, since the difference between 
the two expressions, 



E 

|u| — n 

E 

H=7 



N(v,q) 
L 

N(v,q) 
L 



i 



qvi 



(q + l) 



E 



Ivi 



(q + l) 



(57) 



vanishes due to the fact that V = Yli v l- This com- 
pletes the prove of relationship (|44[) with semiclassical 
means. This result has the added bonus of showing that 
we recover the fidelity amplitude in Eq. (fT]) for r < 1 for 
systems with broken time reversal symmetry as we know 
that the semiclassical and RMT parametric form factor- 
exactly agree in that regime [2ll . |22j . 

To show semiclassically that Eq. (|44l) also holds in the 
symplectic case (f3 = 4) we consider spin-orbit interaction 
for spin 1/2-particles. In this case the spectral form fac- 
tor K(t; 7) on the left hand side of Eq. (|4"4")l is modified to 
-\Kgoe(-t/2;t) with K GO e(t;j) being the form 
factor in the GOE-case. The right hand side yields for 
spin-orbit interaction the additional factor (— 1/2) L ~ V 
that was derived in [32j from the results obtained in [37j • 
A short calculation then shows that Eq. (|4"4"|) also holds 
in the symplectic case. However, we note here that it is 
not true for general spin s-particles. 



IV. EHRENFEST TIME DEPENDENCE OF 
THE FIDELITY AMPLITUDE 

The Ehrenfest time te (see e.g. Ref. (38|) is the 
timescale a minimal wave-packet needs to spread in the 
phase space of a chaotic system to a size such that it 
can no longer be described by a single classical trajec- 
tory. As was pointed out by Aleiner and Larkin [39j in 
the context of transport problem, te is the minimal time 
required for quantum effects to appear. Because of this 
Ehrcnfest-time effects on stationary transport processes 
have been a subject of considerable interest, both theo- 
retically m, US 111 HI] and experimentally [4Q|. Further- 
more, signatures of the Ehrenfest time were also studied 
in the time domain [lj| [3l], H2, 5H |42j , where they are 




FIG. 4: Sketch of trajectories 7 and 7' with a 2-encounter il- 
lustrating the length restrictions in the case of finite Ehrenfest 
time. 



expected to be particularly pronounced. In this section 
we investigate T^-effects on the fidelity decay. Namely, 
we will consider the ^-dependence of the first quantum 
correction to the fidelity amplitude, Eq. (|34|) . for systems 
with time reversal symmetry Our treatment will follow 
the lines of Refs. 0,^. 



Below we derive Ehrenfest time corrections to 



GOE 



coming from the trajectories of the type (A), see Fig. [1] 
To this end we have to specify the classical constant c 
appearing in Eq. ([2"4"| . Note that the different encounter 
stretches are subject to uncorrelated disorder when their 
spatial distance becomes larger than £. We thus require 
that inside of an encounter two segments of the trajec- 
tory are separated by distances less then £, see Fig. 2) 
Accordingly, the encounter time is defined by 



t 



A \CXb \su\ 



(58) 



where C is the system size. The range of validity of our 
approach is given by relation ()37[) . 

The densities ([2"rH) and |[2T[) should be multiplied by a 
Heaviside function in time ensuring that a contribution 
exists only if the trajectory time t is sufficiently long to 
enable a closed path. Specifically, on the right hand side 
of the encounter, depicted in Fig. [4J the stretches should 
be separated by a distance of the order of C in order 
to close themselves and form a loop. On the left side, 
however, the encounter stretches have to be separated 
only by the distance £ (as we consider the case (A)). This 
means the minimal time of the trajectory is 2< onc + 2t^L, 
where 



t $L = A" 1 ln(£/0 



(59) 



is the time it takes the stretches to be separated by the 
distance C when they are initially separated by the dis- 
tance £. 

Accordingly, the weight function ([2"6"| is slightly mod- 
ified by introducing these minimal times and takes the 
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form 



vh,a(u, s) 



(t-2(t cnc + t eL ))' 



9(t-2t enc -2tfr) . (60) 



To account for the correction due to the action difference 
in Eq. (f!?2")) we again use i cnc as defined in Eq. ([55]) . Now 
we are in the position to calculate the Ehrenfest time 
corrections to mQ OE (t; V) in the case (A). As before, we 
will consider the Laplace transform (|45| of the fidelity 
amplitude. With the above corrections it reads 

n (pr) = dte-^^^I(a,t) , (61) 

where we shifted the integration range by 2t onc + 2t^, 
and the integral 



I{a,t) = — 



du i ds- exp (-2at cnc - 4rt cnc ) 

JO ^enc 

(62) 

can be evaluated as in App. B of Ref. [321 ] . In a simi- 
lar way the calculations can be performed for trajectory 
pairs of type B with encounters at the beginning or at 
the end of the trajectory. 

Summing up all the contributions and taking the in- 
verse Laplace transform we obtain the final result for the 
Ehrenfest time dependence of the entire, leading quan- 
tum correction in the orthogonal case: 



m£(t;r) 



(f - 2T d E ) 2 e- T[t+2T ^9{t - 2r|), (63) 



where r% = A^ln^/As), r| = A- 1 ln(^ 2 /(£As)) and 
2t e = t e + t e . Equation (|63| reveals the role of finite 
Ehrenfest times corrections: they lead to a time shift 
such that the quantum corrections set in later, i.e., for 
short times the system behaves "classically" . 

Note also that corresponding calculations for the spec- 
tral parametric correlations give a time delay 2t e which 
is independent of £, since in that case only periodic tra- 
jectories are considered. This demonstrates that in gen- 
eral the connection @ between the fidelity amplitude 
and spectral parametric correlations breaks down when 
non-universal effects such as the Ehrenfest time correc- 
tions are taken in account. 



studies of periodically kicked systems which can be re- 
garded as time-dependent, in Ref. [43[ in the context of 
decoherence and in Ref. [46j . which represents a direct 
extension of the semiclassical approach of Refs. 0, to 
a spatially and time-varying random potential. 

In Ref. [46| , both a finite disorder correlation length £ 
(and associated time ^/vq with vq the particle velocity) as 
well as a finite correlation time To characterizing tempo- 
ral fluctuations was introduced. It was shown semiclas- 
sically on the level of the diagonal approximation, how 
the FGR decay is governed by a decay rate T into which 
both time scales generally enter and which is predomi- 
nantly determined by the shorther of the two times, £/«o 
or To, if they strongly differ. 

In this section we will use the perturbation model of 
Ref. [46| and illustrate the effect of a spatially and time- 
dependent perturbation for the representative case of the 
first quantum correction, m q Q OE (t] T), see Eq. (|34[) for the 
static case. Since this interference contribution is based 
on the mechanism that the same (static) perturbation 
exists along the two encounter stretches traversed at dif- 
ferent times we can compute how finite correlation times 
To will reduce this effect. 

We will consider the interesting case where £/t>o -C To 
so that r depends only on the spatial fluctuations and 
not on To. (For pure temporal fluctuations or £/i>o 3> To, 
a change in To will alter the exponential decay rate T and 
thereby mask the effect of To on quantum fidelity contri- 
butions.) We assume that the time dependence and the 
spatial dependence of the perturbation can be separated, 



i.e. 



AS* 7 = / dt'V(q(t'))W(t'), 



(64) 



which allows for further analytical treatment. We first 
consider the averaged phase difference for this perturba- 
tion 



exp [ — A.S'.y 



exp 



2 



dtdt'V (q(t))W(t) 



XV (q(tj + 2< cnc - t')) W(ti + 2t cnc - t') 



V. FIDELITY AMPLITUDE FOR A 
TIME-DEPENDENT PERTURBATION 

Many physical situations require a generalization of 
fidelity decay to systems with a time-varying perturba- 
tion. This is relevant if a subsystem evolves under the 
influence of a time-dependent external environment or 
may arise if one models fidelity decay for a many-body 
system in terms of single-particle dynamics exposed to 
an external fluctuating potential mimicking the mutual 
interactions. Fidelity decay for time-dependent pertur- 
bations has been addressed in Refs. [44l, foal in numerical 



x exp (— Yt) 



(65) 



with a loop of length ti connecting the two encounter 
stretches. The contribution in the last line of this equa- 
tion is the overall exponential decay, not taking into ac- 
count the correlated way the perturbation acts during the 
encounter. The additional effects of this correlation are 
included in the exponential in the second and third line 
in the equation above, and yielded a contribution pro- 
portional to e _2rt ° nc for an explicitly time-independent 
perturbation. However we now have to analyze this con- 
tribution in more detail: First we again use that the two 
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stretches during the encounter are very close together im- 
plying q(t; + 2t cnc — if) w q(tf). Furthermore, we assume 
in this section, as in Ref. [46|, a Gaussian form of the 
spatial and the time-dependent perturbation 

(V (q(t)) W(t)V (q(U + 2t cnc -t')) W(U + 2t cnc -t')} 



Tv h 2 

17^ 



exp 



|q(t)-q(Or 
£ 2 



x exp 



(tj + 2i cnc — if — ty 



(66) 



The two time integrals in Eq. (|65|) are transformed into 
one integral with respect to r = t — if and one with 
respect tor' = (t+t')/2. The integral with respect to r is 
performed from — oo to oo assuming that the correlation 
length of the spatial part, V (q(t)), is much shorter than 
the one of the time dependent part W(t), £ <C v t . The 
integral with respect to r' is from to i onc , and we thus 
obtain from Eq. (f6l)]) 

(exp (|A5 7 



exp 




x exp 



We perform the two time integrals taking into account 
that only terms linear in t cnc will give a contribution 
when performing the s, it-integrals, see the remark after 
Eq. (f2"9"| . This finally yields for the action difference due 
to the perturbation 



exp 



l-2rcxp U (t O nc+0(tenc)) 



(68) 



exp (-rt) . 



The first term in the squared bracket gives a contribution 
independent of the perturbation and is canceled by the 
contribution coming from case (B). We insert the remain- 
ing term into our expression for mg 0£ (f; I\ To), evaluated 
as m Sec. HH but that now depends additionally on the 
correlation time To of the explicitly time dependent part 
of the potential: 

rt-2t e 



dsdu 



dt. 



(t - 2t c 



x (-2r) e« su exp 



-rt- -V 

T 6 



(69) 

After performing the remaining integrals we finally ob- 
tain the leading quantum correction in the presence of 
an explicitly time dependent perturbation 



2r 



-T iErf 



1 — exp — j 



exp (-rt) 



(70) 



with the error function Erf(x) = -J= f Q dt exp (— t 2 ). By 
expanding Erf(x) we obtain in the limit t <C tq: 



rt 2 

l GOi?(t < T o; r ) - ~^~ e 



-r/ 



1 / t 



, (71) 



indicating a small reduction of the interference term (|34jl , 
— (rt 2 /Tn) exp(— rt) of the static case. 

Much more interesting and relevant is the opposite 
limit, t 3> tq, where we find 



l GOE 



Ft 2 

(t>T ;r)~ -— e 



-r< 



t V * 



(72) 



Quantum fidelity corrections ([34]) in the static case arise 
at time scales t given by the geometrical mean, (t^Tu) 1 / 2 , 
of the Heisenberg time and the decay time tj = 1/21", 
Eq. (fT4"|) . Equation (|T2"|) hence implies, in view of the hier- 
archy of timescales, £/v <C 2t^ < Th, that the quantum 
fidelity contributions are suppressed (by To/t) compared 
to the static case, if r is much smaller than the above 
mentioned geometrical mean that usually represents a 
large time scale. Together with the initial assumption, 
£/i>o <C To, we can conclude that quantum fidelity con- 
tributions are suppressed for time- varying perturbations 
with £/vq <C t -C (^Ih) 1 / 2 . Furthermore, such a sup- 
pression of this negative quantum correction implies that 
upon reducing To, that is introducing faster time varia- 
tions, the overall fidelity amplitude increases in the FGR 
regime and approaches exp(— Ft). 



VI. FIDELITY 

In this section we study the effect of the loop correc- 
tions on the semiclassical expression of the fidelity itself, 
i.e. the average of the squared modulus of the fidelity 
amplitude: M sc (t) — (|m sc (t)| 2 ). We consider both the 
FGR and Lyapunov regime. 



A. Fermi-golden-rule regime 

In Ref. [ll[ the fidelity has been addressed in the RMT 
approach within the linear response approximation valid 
in the transition region between the perturbative and 
FGR regime. Our semiclassical approach is not limited 
to such weak perturbations. ^From Eq. ([7]) we obtain for 
the fidelity the semiclassical expression 
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M sc (t) 



\2irh) 



«/r./r„./q</r./r„./q«>; i ( r + ~) #o (r - |) *o (r + ~J *S f r - | 



EV- ^1/2^1*1/2^*1/2^1/2 
2^ °7l °72 ° 7 £ 

7l(ro->r,t), 72( r 0^ r <*). 



(73) 



x exp 



'72 



■pf 



A5 7l - AS 7s 



As the fidelity decay in the FGR regime arises predom- 
inantly from uncorrelated trajectories 71 and 72 we can 
perform an disorder average of the phase differences AS yi 
and AS J2 independently and obtain an expression con- 
taining four propagators that resembles the one obtained 
for the variance of the survival probability in Ref. [32| . 
with the only difference that the openness of the system is 
replaced by the presence of the random perturbation, i.e. 
Eq. by Eq. (JUJ). In Ref. [13] it could be shown that 
the variance vanishes in the semiclassical limit h —* 0. 
This argument still holds true here. Hence we conclude 
that there remain also only negligible corrections to the 
fidelity apart from those arising from the squared modu- 
lus of the average fidelity amplitude already calculated, 
implying M sc (t) = (\m sc (t)\ 2 ) sa (|m sc (*)|) . Hence we 
can simply use the right hand side of the last relation 
to obtain, to leading order in t/Tn, the following for the 
GOE case 



MqoeW = exp(-2H) ( 1 



and for the GUE case 



2Tf 3 

Mg UE (t) = exp(-2rt)(l--^ 



(74) 



(75) 



B. Lyapunov regime 

Above we considered the case where the disorder av- 
erage in the calculation of M sc (t) can be performed in- 
dependently for the trajectory pairs occurring in the cal- 
culation of each m sc (t). However, already on the level of 
the diagonal approximation there is a further contribu- 
tion originating from configurations where all four tra- 
jectories are too close together to perform the disorder 
average independently. The regime of large T where the 
FGR terms (|74I75|) are rapidly decaying and this contri- 
bution becomes important is referred to as the Lyapunov 
regime, because it decays as e _At [1]. 

Here we examine whether additional contributions may 
arise from loop diagrams. To this end we briefly review 
the semiclassical calculation of the diagonal contribution 
to the fidelity in the Lyapunov regime [4J and consider 
afterwards the role of trajectories differing at encounters. 



Starting from Eq. (|73p . one performs the disorder av- 
erage along the trajectories no longer independently for 
both action differences A5* 7l and A5 72 . For two nearby 
trajectories 71,72, one instead linearizes the motion of 
one trajectory around the other to obtain 



AS 7l - A5 72 = [ dt'VLl (f ) [q 7l (f ) - q 72 (f )] , 
Jo 

(76) 

where L 7l (£') is the Lagrangian associated with the dis- 
order as in Eq. ((SJ), and q 7 . (t') denotes the coordi- 
nates of the trajectories 7$ at time t' . The difference 
q 7 i (*0 — q 72 (*') can t nen be expressed by the difference 
of the final points r and r of 71 and 72, respectively, 
using the exponential separation of long neighboring tra- 
jectories in the chaotic case. Assuming again a Gaussian 
distribution of the random variable, AS 7l — A5 72 , as in 
the FGR regime for A5 7 in Eq. (fTT))) . one obtains 



exp 



-(as 71 -as 72 ; 



(77) 



exp 



.-Ly t ^*dtV( t '+*"- 2 *)a vs (n - r 2 f 



with the force correlator C V s = (V£% {t')VL^ (*")}. 

Replacing in Eq. (|TT[) the second time integral by 
one with respect to t' — t" and taking into account the 
short range behavior of the force correlator (it vanishes 
on scales larger than the correlation length £) one can 
perform the (tf — i")-integral in the range from —00 to 
00. Further evaluation of the i'-integral, where one ne- 
glects contributions from the lower limit because they are 
damped by a factor e _At , eventually gives 



exp 



AS. 



12 , 



exp 



A 



ri - r 2 



(78) 

Here the constant A depends on £ and the disorder 
strength. Equation (|78p is afterwards inserted into the 
full expression for the fidelity finally yielding a contri- 
bution proportional to e _At within the diagonal approxi- 
mation (ij. We note that when performing the two time 
integrals in Eq. (|77|) . only contributions from t' and t" 
close to t mattered. 

We now consider a possible effect of non-diagonal loop 
contributions. The basic contribution of this kind would 
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originate from trajectories 71 and j[ forming pairs as 
depicted in Fig. 1 and two nearly identical trajectories 
72 and 7 2 . When calculating CVs, correlations between 
points of 71 traversed at different times could get impor- 
tant as the latter quantity depends on the difference of 
the positions along 71 at two different times. However, 
as we noted below Eq. (|78p , correlations in the Lyapunov 
regime only matter between points at the end of the orbit, 
i.e. between points, where the motion can still be approx- 
imated by a free motion. Compared to that distance, an 
encounter stretch is exceedingly long, so that correlations 
between different encounter stretches cannot play a role 
for this contribution. This implies that there is no signif- 
icant effect of such orbit pairs in the Lyapunov regime, 
because there exists, apart from the action difference and 
the weight function, no further phase change induced by 
the encounter (depending on the s, u-coordinates). 

The same reasoning can be directly carried over to cor- 
related orbits differing by an arbitrary number of encoun- 
ters, and hence no off-diagonal interference contributions 
from orbits differing in encounters can be obtained in 
the Lyapunov regime. This constitutes another distinct 
difference between FGR- and Lyapunov decay. 

VII. CONCLUSIONS AND OUTLOOK 

In this work, we have presented a semiclassical theory 
for the time decay of the quantum fidelity amplitude and 
the fidelity, or Loschmidt echo, for quantum systems that 
are chaotic in the classical limit. Our focus has been on 
the calculation of quantum contributions to the fidelity 
amplitude beyond the leading term, ~ e~ rt , which we 
showed to arise semiclassically from interference between 
pairs of topologically related trajectories involving classi- 
cal correlations between paths. Applying advanced semi- 
classical methods to deal with multiple pairs of correlated 
orbits, we derived quantum corrections which agree with 
supersymmetric results in the universal Fermi golden rule 
regime. Deriving recursion relations for the semiclassical 
objects, we could show that the semiclassical formulation 
of the fidelity amplitude obeys unitarity and we could 
confirm the interesting relation |(3J) between the fidelity 
amplitude and the parametric form factor to any order 
in t/Tn for times smaller than the Heisenberg time Th- 

Besides, our semiclassical approach provides insight 
into the interference mechanisms underlying quantum fi- 
delity decay. The leading quantum corrections in a t/Tn- 
expansion, Eq. (|34|l . can be regarded as arising from a 
weak-localization type effect leading to a reduction in the 
fidelity amplitude in the time-reversal case, which is sus- 
ceptible to time-reversal symmetry breaking, e.g. through 
a magnetic field. Moreover, the present approach en- 
ables an interpretation of the fidelity dependence on the 
characteristics of the assumed random perturbation, in 
particular the spatial correlation length £ and the time 



correlation t$. The first appears as a further length 
scale in the Ehrenfest time dependence of the fidelity 
amplitude. Even more interestingly, through To we can 
estimate effects from time- or frequency-dependent per- 
turbing sources, such as phonons in solids. As a result, 
we could derive explicit expressions, see Eqs. (|7QI71l72p . 
showing how the negative weak localization-type contri- 
bution is suppressed with decreasing To, and thereby the 
fidelity is counter-intuitively enhanced. 

One main open question arising from this work is how 
to extend the present semiclassical formalism which is 
limited to time scales below the Heisenberg time, to 
longer times. This would also open up a way towards 
a semiclassical understanding of the fidelity amplitude 
revival close to Th in the FGR regime [H, [Hj]. More- 
over, semiclassics beyond Th would allow for extending 
our analysis to a further decay regime, namely for treat- 
ing the perturbative regime of very weak perturbations 
corresponding to energy scales below the mean level spac- 
ing. 

Throughout this work we considered fidelity decay due 
to a global perturbation affecting the whole configura- 
tion space or system boundary. Recently, the comple- 
mentary case of local, not necessarily weak, perturba- 
tions has been assessed semiclassically on the level of the 
diagonal approximation [47], |48j]. There, instead of the 
FGR- and Lyapunov decay, corresponding decay regimes, 
namely a FGR-type and a new 'escape rate' regime, 
have been identified exhibiting in particular a peculiar 
non-monotonic crossover between the two regimes when 
passing from weak to strong local perturbation. Quan- 
tum corrections in both regimes would originate from the 
same diagrams as for semiclassical decay of open quan- 
tum systems, and a semiclassical treatment of the latter 
can be readily applied to obtain a more complete 
picture of fidelity decay for local perturbations. Again, 
accessing time scales beyond the Heisenberg time remains 
there as a challenge for the future. 

Finally, our semiclassical techniques may be applied to 
extend a recent semiclassical analysis [49j of the mech- 
anism of time-reversal mirroring [50], a concept related 
to the fidelity, by including quantum interference effects 
from correlated mirror trajectories. 

VIII. ACKNOWLEDGMENTS 

We thank Inang Adagideli, Peter Braun, Arseni Gous- 
sev, Philippe Jacquod, Rodolfo Jalabert, Heiner Kohlcr, 
Cyril Petitjean, Thomas Seligman and Hans-Jiirgen 
Stockmann for useful discussions and Arseni Goussev, 
Rodolfo Jalabert for a critical reading of the manuscript. 
Financial support by the Deutsche Forschungsgemein- 
schaft within SFB/TR 12 (BG), GRK 638 (DW, MG, 
JK, KR), FOR 760 (KR) and from the Alexander von 
Humboldt Foundation (JK) is gratefully acknowledged. 



14 



[1] A. Peres, Phys. Rev. A 30, 1610 (1984). 

[2] G. Usaj, H.M. Pastawski, and P.R. Levstein, Mol. Phys. 
95, 1229 (1998). 

[3] A. Abragam, The principles of Nuclear Magnetism (Uni- 
versity Press Oxford, 1961). 

[4] H.M. Pastawski and R.A. Jalabert, Phys. Rev. Lett. 86, 
2490 (2001). 

[5] Ph. Jacquod, P.G. Silvestrov, and C.W.J. Beenakker, 

Phys. Rev. E 64, 055203(R) (2001). 
[6] N.R. Cerruti and S. Tomsovic, Phys. Rev. Lett. 88, 

054103 (2002). 

[7] T. Prosen, T.H. Seligman, and M. Znidaric, Prog. Theor. 

Phys. Suppl. 150, 200 (2003). 
[8] F.M. Cucchietti, H.M. Pastawski, and R.A. Jalabert, 

Phys. Rev. B 70, 035311 (2004). 
[9] T. Gorin, T. Prosen, T.H. Seligman, and M. Znidaric, 

Phys. Rep. 435, 33 (2006). 
[10] For a review of related quantities describing e.g. decoher- 

ence, see P. Jacquod and C. Petitjean, Adv. in Physics 

58, 67 (2009). 

[11] T. Gorin, T. Prosen, and T. Seligman, New. J. Phys. 6, 
20 (2004). 

[12] H.-J. Stockmann and R. Schafer, New J. Phys. 6, 199 
(2004). 

[13] H.-J. Stockmann and R. Schafer, Phys. Rev. Lett. 94, 
244101 (2005). 

[14] H. Kohler, I. Smolyarenko, C. Pineda, T. Guhr, F. 
Leyvraz, and T.H. Seligman, Phys. Rev. Lett. 100, 
190404 (2008). 

[15] J. Vanicek and E. Heller, Phys. Rev. E 68, 056208 (2003). 
[16] W. Wang, G. Casati, B. Li, and T. Prosen, Phys. Rev. E 

71, 037202 (2005). 

[17] M. Sieber and K. Richter, Physica Scripta T90, 128 

(2001) . 

[18] S. Miiller, S. Heusler, P. Braun, F. Haake, and A. Alt- 
land, Phys. Rev. Lett. 93, 014103 (2004); Phys. Rev. E 

72, 046207 (2005). 

[19] P.W. Brouwer, S. Rahav, and C. Tian, Phys. Rev. E 74, 
066208 (2006). 

[20] S. Heusler, S. Miiller, A. Altland, P. Braun, and F. 

Haake, Phys. Rev. Lett. 98, 044103 (2007). 
[21] T. Nagao, P. Braun, S. Miiller, K. Saito, S. Heusler, and 

F. Haake, J. Phys. A 40, 47 (2007). 
[22] J. Kuipers and M. Sieber, J. Phys. A 40, 935 (2007). 
[23] D. Waltner, S. Heusler, J.-D. Urbina, and K. Richter, J. 

Phys. A 42, 292001 (2009). 
[24] K. Richter and M. Sieber, Phys. Rev. Lett. 89, 206801 

(2002) . 

[25] i. Adagideli, Phys. Rev. B 68, 233308 (2003). 

[26] S. Heusler, S. Miiller, P. Braun, and F. Haake, Phys. Rev. 



Lett. 96, 066804 (2006). 
[27] P.W. Brouwer and S. Rahav, Phys. Rev. B 74, 075322 
(2006). 

[28] Ph. Jacquod and R.S. Whitney, Phys. Rev. B 73, 195115 

(2006) . 

[29] G. Berkolaiko, J.M. Harrison and M. Novaes, J. Phys. A 

41, 365102 (2008). 
[30] for an introduction into these semiclassical techniques 

see, e.g.: D. Waltner and K. Richter, in Nonlinear 

Dynamics of Nanosystems, edited by G. Radons, H.G. 

Schuster, and B. Rumpf (Wiley- VCH, Berlin, 2010). 
[31] D. Waltner, M. Gutierrez, A. Goussev, and K. Richter, 

Phys. Rev. Lett. 101, 174101 (2008). 
[32] M. Gutierrez, D. Waltner, J. Kuipers, and K. Richter, 

Phys. Rev. E 79, 046212 (2009). 
[33] J. Kuipers, D. Waltner, M. Gutierrez, and K. Richter, 

Nonlinearity 22, 1945 (2009). 
[34] K. Richter, D. Ullmo, and R.A. Jalabert, J. Math. Phys. 

37, 5087 (1996). 
[35] M. Gutzwiller, Chaos in Classical and Quantum Mechan- 
ics (Springer Verlag, New York, 1990). 
[36] M. Gutierrez and A. Goussev, Phys. Rev. E 79, 046211 

(2009). 

[37] J. Bolte and D. Waltner, Phys. Rev. B 76, 075330 (2007). 
[38] B.V. Chirikov, F.M. Izrailev, and D.L. Shepelyansky, 

Sov. Sci. Rev. C 2, 209 (1981). 
[39] I.L. Aleiner and A.I. Larkin, Phys. Rev. B 54, 14423 

(1996). 

[40] O. Yevtushenko, G. Liitjering, D. Weiss, and K. Richter, 

Phys. Rev. Lett. 84, 542 (2000). 
[41] H. Schomerus and J. Tworzydlo, Phys. Rev. Lett. 93, 

154102 (2004). 

[42] H. Schomerus and P. Jacquod, J. Phys. A 38, 10663 
(2005). 

[43] F.M. Cucchietti, D.A.R. Dalvit, J.P. Paz, and W.H. 

Zurek, Phys. Rev. Lett. 91, 210403 (2003). 
[44] Ph. Jacquod, P.G. Silvestrov, and C.W.J. Beenakker, 

Phys. Rev. E 64, 055203(R) (2001). 
[45] G. Benenti and G. Casati, Phys. Rev. E 65, 066205 

(2002). 

[46] F.M. Cucchietti, C.H. Lewenkopf, and H.M. Pastawski, 

Phys. Rev. E 74, 026207 (2006). 
[47] A. Goussev and K. Richter, Phys. Rev. E, 75, 015201(R) 

(2007) . 

[48] A. Goussev, D. Waltner, K. Richter and R.A. Jalabert, 

New J. Phys. 10, 093010 (2008). 
[49] H.L. Calvo, R.A. Jalabert, and H.M. Pastawski, Phys. 

Rev. Lett. 101, 240403 (2008). 
[50] C. Draeger and M. Fink, Phys. Rev. Lett. 79, 407 (1997). 



